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Abstract 

Traffic flow at low densities (free traffic) is characterized by a quasi-one-dimensional 
relation between traffic flow and vehicle density, while no such fundamental diagram 
exists for 'synchronized' congested traffic flow. Instead, a two-dimensional area of 
widely scattered flow-density data is observed as a consequence of a complex traffic 
dynamics. For an explanation of this phenomenon and transitions between the dif- 
ferent traffic phases, we propose a new class of molecular-dynamics-like, microscopic 
traffic models based on times to collisions and discuss the properties by means of 
analytical arguments. Similar models may help to understand the laminar and tur- 
bulent phases in the dynamics of stock markets as well as the transitions among 
them. 
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1 Introduction 

The dynamics of vehicle traffic has challenged researchers for more than five 
decades now. Some major breakthroughs have been made only recently thanks 
to physicists who applied methods from statistical physics and non-linear dy- 
namics. In the meantime, empirical observations indicate that there are dif- 
ferent phases of traffic. According to Kerner [1-8], these are (see Fig. 1): 

1. free traffic (FT) characterized by a unique flow-density relation Q e (p) up to 
some maximum flow Q max , 



Preprint submitted to Elsevier Preprint 



1 February 2008 



2. wide moving clusters (traffic jams) characterized by a jam line J(p) and 
'natural constants' such as the propagation velocity C < and the outflow 
Qout < Qmax °f wide clusters, and 

3. 'synchronized' congested flow, characterized by a synchronization of average 
vehicle speeds among neighboring lanes and a wide scattering of flow-density 
data. 

Synchronized flow has three subtypes: (i) stationary and homogeneous states, 
(ii) oscillatory states, and (iii) homogeneous-in-speed states, where the velocity 
is constant, i.e. the flow is proportional to the density. While the transitions 
among these subtypes are continuous, the transitions between free traffic, wide 
moving jams, and synchronized flow are hysteretic in nature [4-8]. 

FlowQ 



Qout 




Pout Per Density p 



Fig. 1. Schematic illustration of measurements of flow-density data (one-minute 
averages) after Kerner. 

Most simulation models show different traffic states as well [9-14]. Usually, 
one finds free and stable traffic at low densities. At higher densities, traffic 
is unstable in most models, which can give rise to emergent traffic jams due 
the slow relaxation to some density-dependent equilibrium velocity V e (p). The 
delayed adaptation leads to overbraking and chain reactions of the following 
drivers, finally resulting in a standstill of cars. In addition to traffic jams, these 
models can reproduce congested traffic states of homogeneous and oscillatory 
type behind bottlenecks reminding of 'synchronized' congested traffic of types 
(i) and (ii) [15,16]. The resulting state depends on the flow on the freeway and 
on the bottleneck strength, but it is also history-dependent, as traffic flow can 
be multi-stable. Despite of this success, some people believe that the obser- 
vation of the wide scattering of congested flow-density data has not yet been 
fully understood. Many researchers in the field of traffic physics have, there- 
fore, published their own theories of 'synchronized' flow, but none of them 
is generally accepted, and the on-going debate on this hot topic is still very 
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controversial. The proposed explanations include mixtures of different vehicle 
types (cars and trucks) [17], a heterogeneity in the time headways [18], changes 
in the behavior of "frustrated" drivers [19], anticipation effects [20-22], non- 
unique equilibrium solutions [23,24], and multiple metastable oscillating states 
[25] . Only the hypotheses of mixed vehicle types and heterogeneous time head- 
ways have been empirically supported. Here, we will suggest another approach 
to the subject, which is empirically motivated. Although we presently do not 
share some of Kerner's interpretations of empirical observations and his criti- 
cism of the existing traffic models, this paper aims at developing a new class 
of models consistent with his interpretation of observed traffic data. 



2 The Times-to-Collision Model 

In the following, we will restrict to identical driver-vehicle units, but distin- 
guish between two different driver behaviors corresponding to different traffic 
states: 

1. In free traffic, a vehicle % located at place Xi(t) at time t will try to adapt 
its actual velocity i>j = dxi/dt to the so-called optimal velocity v e (Axi) = 
V e (l/Axi) within a certain relexation time: 



v e (Ax) is a function of the headway (gross/brutto distance) Axi = (xi-i — 
Xi) to the vehicle in front and can be empirically determined. Alternatively, 
we may use the simple relation 



where Vq is the desired (maximum) velocity in free traffic, I is the minimum 
vehicle distance (effective vehicle length) and T the safe time clearance 
(net/netto time gap). This corresponds to the idea that a vehicle tries to 
keep a velocity-dependent safe distance (/ + Tv) and to drive at maximum 
speed vq when there is enough free space in front. 
2. In 'synchronized' congested traffic, vehicles cannot move freely anymore. 
Their main concern is to avoid accidents. That is, the time to collision 
(TTC) given by Tj = (Axj — l)/Avi, where Avi = (i>j — v j_i) is the relative 
velocity, should not drop below a certain desired value T , in accordance with 
observed driver behavior. In other words: Vehicles tolerate a small (netto) 
clearance (Axj — I) to the vehicle in front, if the relative velocity Av^ is 
small. Moreover, the headway to the leading vehicle is history-dependent. 
It may suddenly shrink due to a lane-changing vehicle, but usually this does 
not cause panic braking in order to restore the safe distance (I + Tv). 



dvj _ v c (Axi) - Vj 
dt r 





(2) 
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Let us first discuss the implications of the empirically justified observation that 
drivers try to maintain a constant inverse time to collision l/Tj = Ai>j/ (Arcj— I) 
during their following behavior. This implies the equation d(l/Ti)/dt = 0, 
which finally results in 

dvj _ dvj_i _ Avj 2 

dt ~ dt C Axi -I ^ ' 



with c = Q(Avi). The Heaviside step function Q(Avi) is 1 when Avi > 
and otherwise 0, thereby reflecting that a vehicle would not brake, when the 
relative velocity Avi is zero or negative (given this is not recommended by a 
deceleration dv^/dt < by the leader). 

We will find that, similar to some other traffic models [10,26-30], the Avr 
dependent term guarantees safe driving. Introducing the new variable Zi = 
(Axi — l) and considering Avi = —dAxi/dt = —dzi/dt, we obtain the equation 

d 2 Zi dAvi c { dzi\ 2 Avidzi 
dt 2 dt Zi \dt J Zi dt 



or 



dhi(Avi/v*) dAvi/dt dzi/dt dhi(zi/z*) 



dt 



Av,: 



dt 



(5) 



which is analytically solved by 



A Vi (t) 
Avi(0) 



Zjjt) 

Zi(0) 



(6) 



We can immediately see that for any c > 0, the relative velocity Avi will 
become zero when Zi = 0, i.e. accidents can be avoided for sure. For c = 1, we 
find 



dAvi 
dt 



Aui(0) 



Axi(0)-l 



Avi, 



(7) 



while the case of kinematic braking with constant deceleration 



dAvi 
dt 



= -b = 



[A«i(0)] s 



2[A^(0) -I] 



< 



(8) 



corresponds to c = 1/2, as usual. Moreover, the equilibrium solution of equa- 
tion (3) is Avi = 0, but the distances Axi between the vehicles in that 
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state are history-dependent and can be very different. That is, maintaining 
a certain time to collision results in driving at some constant velocity as in 
homogeneous-in-speed states, but the local density p = 1/Axj may vary con- 
siderably. 

For c = 1, equation (4) is solved by 

z i (t) = z i (0)exp(-^^t) , (9) 



which confirms the history-dependent relaxation behavior. If Avi > 0, the 
headway exponentially approaches the minimal distance Axi = I. Therefore, 
we expect platoon formation of fast cars behind slow ones with Avi < 0. In 
contrast, the headways in front of slow vehicles will usually become so high 
that they will eventually switch to free driving (see Sec. 4). Therefore, our 
model assumes driver reactions to the acceleration dvi-i/dt of the leading 
vehicle only during following behavior, when the headways are small. As it is 
hard for drivers to measure accelerations, it would be reasonable to replace 
the acceleration of the leading vehicle by some step function of dvi-i/dt. 



3 Variants of the Model 



We will now introduce an extension of our model of following behavior in 
'synchronized' congested traffic. Let us assume that vehicles slowly tend to 
adapt their time to collision Tj = (Ax, — I)/ Av^ to some preferred value T . 
We may describe this by the equation 



dvi dvi-i Avi 
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-c . 1 1 +u\Ax l -l-T Av i ) (10) 



dt dt Axi — I 
with a small "oscillation frequency" uj. This equation is equivalent to 

+ u z T — + u z Zi = — — , 11 



dt 2 dt Zj \ dt 



which can be analytically solved. We can distinguish the following cases: 

1. For c ^ 1, we can apply the transformation = y]^ 1- ^-, which results in 
the equation of a damped oscillator: 

^WT f + ( l-c)^0. (12) 
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The general solution is, therefore, of the form 



t (t) = (A ie ^ - A 2 e-^) 1/(1 ~ C) , (13) 



where A\ > A 2 are given by the initial conditions and a± — (a — u 2 T /2), 
a 2 = (lu 2 T /2 + a) with a 2 = [(lu 2 T /2) 2 - 7c; 2 ]. When a 2 is negative, the 
solution describes a damped oscillation, otherwise the relaxation to Zi = 
is overdamped. In the case a — 0, the solution is 



(14) 



with parameters B 1 and B 2 given by the initial conditions. 
2. For c = 1, the solution reads instead 



Zi(t) = Ciexp 



(^ + C 2 e- w2Tot ) 



(15) 



co 2 T 

The parameters Ci and C 2 > are again given by the initial conditions. 

In conclusion, we find damped or overdamped non-linear oscillations. The 
system oscillates around and/or converges to Avi = and = 0, i.e. Axi = I. 
Not only would this imply the possibility of accidents, it also means that 
the resulting state is characterized by one or several dense vehicle platoon(s). 
Instead of this, we would prefer oscillations around the safe distance (/ + Tvi) 
and Avi = 0. Therefore, we modify equation (10) a little: 



dvi dvi-i Av. 



2 



dt dt Axi — I 



-c- — l —+u\Axi-l-Tvi-T Avi). (16) 



Between Eqs. (10) and (16) lies a major difference. While there is no stationary 
solution of equation (10) (Avi always tends to be non-zero corresponding to a 
pushy driver behavior), Eq. (16) has the stationary solution V; L = (Arrj— l)/T. A 
linear stability analysis indicates that, in the presence of small perturbations, 
we can expect damped or overdamped oscillations. Instead, we would prefer 
an unstable solution leading to permanent oscillations in the system due to 
vehicle interactions (see Fig. 2). In that case we could expect to find a wide 
scattering of flow-density data in the congested regime based on a complex 
dynamics. 

We can reach this by generalizing equation (11) according to 

^ + ^ W( ,_ r „, H £(^ 2 (17) 
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Fig. 2. Measured oscillations of the clearance z = (Ax — I) and of the relative velocity 
Av around Av = (cf. Hoefs, 1972) indicate an instability in car-following behavior 
(after Helbing and Tilch, 1998). 

with the a friction function assuming negative values in the neighborhood of 
(Axi -l)= Vi T, e.g. 



fx(zi) = u 2 T + to 







( 



Tvi 



+ 



-7 



(18) 



The term in the square brackets is new and causes a negative friction effect, i.e. 
a driving effect, if fi < 0. Such negative friction effects are rather common in 
models of traffic and self-driven many-particle systems [14]. By construction, 
the driving effect can appear in the neighborhood of Zi = v{T only if (u)Tq+2(3— 
7) < 0. The larger (7 — 2(3 — ujT ) is, the larger will the emerging oscillations 
be. Because of the scaling of Zi by v{T, they will also grow with increasing 
Vi, i.e. with local decreasing density p ~ 1/(1 + Tv/). A similar result can be 
found for related approaches for the friction function or for oo 2 itself. 

Under the above assumptions, the oscillations are limited in size, as the fric- 
tion coefficient becomes positive, when the deviation from Zi = v{F becomes 
significant. Moreover, the oscillations are non-linear and influenced by the rel- 
ative velocity Avi = —dzi/dt as well. The corresponding equation describing 
the following behavior of driver i reads 

dvi dvi-i Avf + r} 2 2 

~dt = ~~dt ° Ax- - I ~ ^ ~ ' U ^ ~~ ~~ ^ ' 



with i] = 0. Setting 77 7^ allows to take into account errors in the estimation 
of Avi or unexpected variations of the relative velocity due to lane changes. 
Moreover, it would certainly be realistic to add some random noise ^(t) to the 
acceleration equations, which will be treated in a forthcoming paper. Another 
reasonable model variant would result by using the acceleration a = u 2 (Axi — 
I) as a model parameter rather than uj. However, the above model shows 



interesting behavior also in the special cases given by c = 1, w = 0, w = 
1/ (tT), P = 7 = 0, or rj = 0. The other parameters are known from empirical 
investigations. 



4 Transitions Between the Different Traffic States 



It is well-known that the model (1) for the driver behavior in free traffic 
produces linearly unstable traffic flow and jam formation when the condition 

dAx 2r v ; 



is fulfilled. We denote the lowest density p = 1/Ax, for which this condition 
is fulfilled, by p cr (critical density). At this density, the flow Q e (p) = pv e (l/p) 
reaches its maximum Q max . 

With some didactically justified simplifications, one can say the following: 
Once jams are formed, they are characterized by a jam line 

J(P) = ^ (21) 



related to v — v e (Ax) = (Ax — l)/T, where C = —l/T corresponds to the 
backward propagation velocity of the traffic jam. The lowest density at which 
the jam line J(p) cuts the function Q e {p) defines the characteristic outflow 
Qout from a traffic jam and the density p out of free flow downstream of it 
(see Fig. 1). Between the densities p out and p cr , traffic is metastable, i.e. a 
perturbation will grow, if its amplitude exceeds a certain critical amplitude 
A(p), otherwise it will fade away. The critical amplitude A(p) decreases with 
growing density and becomes zero at p = p a . At this density, the probability 
of a breakdown of traffic is 100%. 

Transitions from free to 'synchronized' congested flow occur more frequently 
than transitions from free traffic to wide moving clusters. We will assume that 
a driver-vehicle unit switches between free driving and following behavior with 
probability 1, when the local flow Q = Vi/Axi reaches the line 

K ^ Q -^-W=j^' (22) 
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i.e. when the time clearance (Axi — l)/vi becomes smaller than 

T' = . (23) 

W max 



The transition probability shall be zero below the jam line J(p), i.e. when the 
time clearance is larger than T. Between the lines J(p) and K(p), we may 
assume that the transition probability P is approximately given by 



with some parameter e > 1. This transition is of hystertic nature. 

Finally, we mention that we have checked out deterministic transition rules as 
well, e.g. a switching from free to following behavior for Axi < (1+T'vi+TyAvi) 
and back to free driving for Axi > (/ + Tvi + T^Avi), where the parameter Tq 
may agree with T . The results will be discussed in another paper. 



5 Summary and Outlook 



We have proposed novel microscopic traffic models assuming that the fol- 
lowing behavior in 'synchronized' congested flow is different in nature than 
free driving. The main assumption is that drivers try to maintain a certain 
time to collision to guarantee accident-free driving. The basic model variant 
had no unique solution, i.e. no fundamental diagram, and the finally evolv- 
ing traffic flow was history-dependent. An extended model variant had a sta- 
tionary solution, but an unstable one, if the parameters were appropriately 
specified. In that case, the vehicle velocities were non-linearly oscillating with 
medium-sized amplitudes around the stationary solution. The conditions for 
transitions between free driving and following behavior were specified in a way 
inspired by the hypotheses on congested traffic by Kerner [4-8]. In particular, 
we should have hysteretic transitions between free traffic, wide moving clus- 
ters, and 'synchronized' flow. Transitions to traffic jams should be rare and 
require a triggering by sufficiently large perturbations. 

Based on our present results, all of this appears to be reproducible with 
the proposed models. Corresponding simulation results will be presented in 
a forth-coming paper. Free traffic is characterized by a one-dimensional flow- 
density relation of positive slope, wide moving clusters by a linear jam line of 
negative slope, and 'synchronized' congested flow by a two-dimensional varia- 
tion of flow-density data. 
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We believe that the transition between laminar regimes (free traffic) and 
strongly varying regimes ('synchronized' flow) may have some analogies with 
the dynamics of stock markets, where one observes time periods of low and of 
high volatility [31,32]. The models suggested above may, therefore, also help 
to gain a better understanding of stock markets. 
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